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In this work we explore how nonhnear modes described by a dispersive wave equation (in our 
example, the nonlinear Schrodinger equation) and localized in a few wells of a periodic potential 
can act analogously to a chain of coupled mechanical oscillators. We identify the small-amplitude 
oscillation modes of these 'coupled wave oscillators' and find that they can be extended into the 
large amplitude regime, where some can 'ring' for long times. We also identify prototypical case 
examples of more complex dynamical behaviour that can arise in such systems including the break- 
down of Josephson-like oscillations and of internal modes more generally, the transfer of energy 
out of/destabilization of fundamental oscillation modes and the emergence of chaotic oscillations for 
large amplitude excitations. We provide details of the phase perturbations required for experimental 
observations of such dynamics, and show that the oscillator formalism can be extended to predict 
large amplitude excitations in genuinely two-dimensional configurations. 



Since Huygens first noticed the synchronization of two 
pendula attached to the same wall [l|, the coupling of 
mechanical oscillators has fascinated scientists and the 
general public alike. Complex and often counter-intuitive 
results are found in even simple coupled systems such as 
the spring-pendulum [2] and two coupled pendula [35 . 
With the advent of exotic coherent wave systems the os- 
cillator formalism has found new life, providing a deep 
understanding of effects such as the Josephson tunnelling 
of supercurrents [4], the evanescent wave tunnelling be- 
tween two optical waveguides [5] and Josephson-like os- 
cillations in a Bose-Einstein condensate (BEC) confined 
in a double- well potential [6|, [D]. Often, the oscillatory 
dynamics in these coherent wave systems are captured 
by the dynamics of a single mechanical oscillator y, |8| . 
In contrast, we explore here how matter waves or optical 
waves in a few wells of a periodic potential, such as light 
in coupled opticall y- in duced waveguides [9| , or a BEC in 
an optical lattice [lO|, may behave analogously to cou- 
pled mechanical oscillators. By focusing on a few (more 
than one) degrees of freedom, we predict new macro- 
scopic wave dynamics, and so open up the possibility of 
experimental investigation of coupled effective oscillator 
dynamics in a precisely controlled environment. 

Tunnelling between multiple wells has been observed 
in BEC in a periodic potential [ll|, however our work 
differs by beginning with a stationary configuration of 



occupied sites in a periodic potential [12|, ll3| , and exam- 
ining the tunnelling dynamics which may occur between 
these nonzero ("excited") sites. Rather than trying to 
characterize all the possible dynamics, as has been the 
emphasis in double- well [IJ] and triple- well 15|, ll6| stud- 
ies, we focus on the periodic excitations, analogous to 
the modes of a coupled-oscillator system. This allows 
us to go beyond a triple-well configuration and identify 
the amplitude/phase oscillation modes which may be ex- 
cited by appropriate initial conditions for multi-site exci- 



tations. We also explore a number of instability pathways 
possible in the multi-well configurations and extend the 
oscillator formalism from a one-dimensional (finite) chain 
to a truly two-dimensional arrangement. We note that 
a recently published work [17| complements our results 
with the detailed derivation of a model for the A^-well 
case and an examination of the stationary states in four 
coupled wells. 

As our prototypical case example, we use systems de- 
scribed by the two-dimensional nonlinear Schrodinger 
equation with a repulsive nonlinearity. 



dt dx'^ dy'^ 



F(x,7/)^-|^P^ = 0, (1) 



and a spatially periodic external potential of depth Vq 
given by V{x^y) = Vb(sin^(x) + sm^{y)). This equation 
appears in different physical contexts, with the most rel- 
evant that of a Bose-Einstein condensate (BEC) in an 
optical lattice [l8|(see e.g. p^] for an example of a con- 
version between the normalised equation ([!]) and phys- 
ical parameters in BEC), or a coherent light field in a 
defocusing material with an optically induced periodic 
potential [9[ . We are interested in the possibility of cou- 
pled coherent wave oscillators described by this model, 
and as a basis for this study we examine the truncated 
nonlinear Bloch waves supported by the periodic poten- 
tial [12|. These states exist inside the linear band-gap 
of the periodic potential through a balance between the 
nonlinearity and Bragg refiection, and form distinct fam- 
ilies depending on the number of occupied sites [19]. In 
the two-dimensional potential the configuration of the 
occupied sites can assume complex forms through suit- 
able contours [l3[ , however for the purposes of this work 
we focus primarily on a finite chain of occupied sites. In 
Fig. [TJa) we show an example of a five site state undergo- 
ing the large amplitude oscillatory dynamics of interest 
to this work. With an appropriate phase perturbation 




FIG. 1: (color online), (a) Schematic of a five-site state lo- 
calized in a periodic potential, showing large scale periodic 
excitation; (b) Appearance of regular nonlinear wave oscilla- 
tions following initial phase perturbation. The colorbar on 
right shows density /intensity according to color. Parameters: 
^ = 6, Vb = 6, e = 7r/4. 



applied to the occupied sites the initially almost equal 
amplitude state starts to periodically exchange particles 
in a coherent and long-lived oscillation (shown in Fig. 
[TJb)). While this dynamics is shown for the model ([T]) 
we note that the primary results of our work are inde- 
pendent of the specific physical details of the problem, 
such as the form of the potential or the precise nature of 
the underlying infinite dimensional wave equation. 

To better reveal the generality of our results we be- 
gin by considering the discrete version of model ^ ob- 
tained from the ansatz ^ 



^M 



jiit))(j)i{x^y)^ where (piix^y) captures the spatial depen- 
dence of the wavefunction at site i for a stationary state 
with M occupied sites and chemical potential ji. Substi- 
tution of this ansatz into ([T]) and integration over the spa- 
tial dimensions leads to time-dependent equations for Ni 
and Oi , the particle number and phase respectively of site 
i. See [1^ and [15| for the derivation in the two-site and 
three-site cases respectively; a derivation for the general 
M site case has recently appeared in [17|. Moreover, such 
models can be considered as finite dimensional trunca- 
tions of the infinite chain cases examined e.g. in 20|, |21|. 



effect of the next order nonlinear coupling terms using 
a nonlinear basis, applicable for the strong coupling be- 
tween sites seen in our continuum stationary states. The 
resulting equations can be found from the Hamiltonian, 



M-l 



H=Y, [2(i^.,^+i - m^^+l{N, + iV,+i))v/iV^iV^ (2) 



cos(i9i+i - Oi) - Xi,i+iNiNi^i (2+ 



M 



However, unlike these earlier works, we consider also the 



cos(2(^,+i - 0,)))] - Y, ^Nf + E,N,, 
using Hamilton's equations for the canonical coordinates 

{0^,Ni), 

As the total mass is conserved, "^^ Ni = TV^, we may set 
Nm = ^T — Xli ~ ^i, leaving M — 1 free mass variables. 
Similarly we have only M — 1 free phase variables, defined 
by the phase differences (6>i+i — 6i). The system (J3j) is 
thus equivalent to a mechanical system with M — 1 de- 
grees of freedom, or more precisely a set of M — 1 coupled 
momentum-shortened mechanical pendula with phase 
(^^+1 — Oi) and conjugate momentum Ni^ and higher- 
order nonlinear coupling given by the rj and x terms. 
The system coupling parameters are given by i^i,i+i = 
- X!^ V(/)i V0i+i + (j)iV{r)(j)i^i - ii(j)i(j)ij^idY/Pi, rji^i^i = 

S-oo^l^i+^^^/^i ^^^ X^.^+1 = I-oo^^^hi^^/^i^ where 
Pi = J_ (pfdr. The effective on-site nonlinearity is 
given by 7^ = J_ (pfdr/Pi and the on-site energy, equal 
between occupied sites in our case, is given by Ei = 
/_ (V(/)i)^ + (/)iF(r)(/)^ — jj.(j)'fdr/Pi. Here the subscripts i 
and i-\-l refer to modes localized in adjacent wells i^i -\- 1 
of the periodic potential. Hence, for instance, in order to 
obtain these parameters in the restricted two-mode case, 
we use 01,2 = (^s±^a)/(2v^) where ^5 and ^a are the 
symmetric and anti-symmetric two-site solutions to the 
full system ([I]) and Ns = Jy^ |(^s + ^a)/2pG^r is the on- 
site mass obtained by integrating over a single occupied 
well Vi of the periodic potential V{x^y). For the case of 
the continuum stationary state with Vb = 6 and jj. = 6 
the associated values are K = 0.17337, 77 = 0.00567, 
X = 0.00085, 7 = 0.196, Nq = 7.64 and E = 1.519 (the 
same for both sites). 

The free boundary conditions of the oscillator chain 
suggest that the ground state occupied site configuration 
does not have equal particle numbers at each site (see the 
discussion for the three-site case [15] and also the results 
of [17] for the four site one). However, we are interested 
in the oscillation modes of the ground state rather than 
details of the ground state itself, so we make the approx- 
imation that Ni = No for all sites (exact in the limit 
K ^ 0) and that all parameters are the same as in the 
two-site case. These approximations allows us to obtain 



explicit forms for the mode structures and frequencies in 
the discrete model. To find these modes, and their associ- 
ated frequencies, we substitute Ni = No-\-Si^ Oi = e^ into 
the system equations (jSj and keep only the linear terms 
in (Si^Ci). We find that the stationary state is neutrally 
stable, as is expected for our defocusing nonlinearity; in 
the focusing case all the norm would tend to focus to 
a single well rendering the above state highly unstable 
with M — 1 real eigenvalue pairs; in the defocusing case, 
there are M — 1 imaginary eigenvalue pairs (correspond- 
ing to eigenfrequencies of oscillation). We list the mode 
eigenfrequencies ujj in Table HI where j = 1 ... M — 1 runs 
from the lowest to highest frequency eigenmode. The an- 
alytical expressions are lengthy in the general case, so for 
the most part we present only the case with ^ = X = 0- 
The calculation of the eigenvalues in the case ?7, X 7^ 
amounts to the diagonalization of a 2M x 2M matrix, 
which for 2 to 5 sites is a straightforward task. The 
eigenvalues for the state /i = 6, Fq = 6 are shown as solid 
black circles in Fig. O 



Two-site cji = yipT - 2r]No)jNo + 4(K - 2r]No) 



Three-site 



Four-site 



Five-site 



0J2 = y^6K-fNo + 9K^ 



CJl 



(4 - 2V2)K-fNo + (6 - 4V2)K^ 

y ^4KjNo + Ik^ 

CJ3 = a/(4 + 2V2)K'jNo + (6 + 4V2)K'' 



UJ2 



UJ2 
CJ3 
UJ4. 



(3 - V^)K-fNQ + ^^'l^^ K^ 
(5 - Vb)K-fNQ + 



(15-5^) 



i^2 



(3 + V^)K-fNQ + (I±|^K2 
(5 + V^)K-iNQ + IH^^^i^2 



TABLE I: Approximate mode eigenfrequencies for multi-site 
stationary states in the weak coupling limit. 

We return now to our original continuum model ([T]) 
and perform a similar linear stability analysis to deter- 
mine the eigenfrequencies and eigenmodes. We see in 
Fig. [2] for the cases of M = 2 ... 5 wells, a confirmation 
of the presence of the M — 1 oscillation eigenfrequencies 
(numbered in ascending order of frequency), as found for 
the discrete case. We see that the discrete model pre- 
dictions are excellent in the two site case, and still very 
good with a larger number of sites despite the assumption 
of equal well populations. This shows that even in the 
strongly coupled case the vast simplification of the dis- 
crete model is still meaningful. Included also in Fig. [2] is 
a schematic representation of the associated eigenmode 
profiles found in both the discrete and continuum cases 
(filled circles). The arrows indicate the direction of am- 
plitude change for one half of an oscillation. In the other 
half the directions of the arrows is reversed. The mode 
structure shows interesting parallels to the more famil- 
iar mechanical oscillator normal modes. Noting that the 




FIG. 2: (color online). Internal mode oscillation frequency 
versus chemical potential fi found from linear stability anal- 
ysis for (a) two-site; (b) three-site; (c) four-site; and (d) five- 
site states. Solid black circles show eigenvalues from the dis- 
crete system at /i = 6. Insets show relative change in den- 
sity/power between sites for each internal mode (filled cir- 
cles with arrows indicating direction of change) and analogous 
coupled effective pendulum modes (open circles). 



symmetry of the phase oscillations mirrors the symmetry 
of the amplitude oscillations, and remembering that the 
phase of the equivalent mechanical oscillator is given by 
(^^+1 — Oi)^ we can represent schematically the analogue 
mechanical mode associated with each amplitude/phase 
oscillation mode. These are shown on Fig. [2] as M — 1 
open circles for each mode, with the arrows now indicat- 
ing the direction of oscillation of the effective pendulum. 
It is interesting to point out here that this effective nor- 
mal mode approach may be more broadly applicable. For 
instance, the works of [22| suggest that the dynamics of 
dark solitary waves as longitudinally vibrating particles 
in BECs share similar symmetries and characteristic os- 
cillation modes. 

An important consideration for the identified oscilla- 
tion modes is their behaviour in the presence of non- 
linearity. In the discrete case the momentum-shortening 
and nonlinear coupling of the effective pendula introduces 
a complexity absent in similar mechanical oscillator sys- 
tems. A key question therefore pertains to whether these 
oscillation modes persist in the full nonlinear wave dy- 



namics (and hence beyond the associated hnearization) . 
And if so, how can they be excited and what role do 
they play in the system's nonlinear dynamics ? The full 
analysis of these modes in the nonlinear case is beyond 
the scope of this work. Here, we restrict ourselves to 
some vignettes showcasing the relevance and potential 
observability of such modes and the intriguing nonlin- 
ear dynamics that they may exhibit. Our approach to 
generation is to look for mode symmetries which are pre- 
served in the presence of the full system nonlinearity. For 
instance one obvious symmetry is a reflection about the 
mid-point of the oscillator chain. Mode 2 in the three-site 
case, mode 2 in the four-site case and modes 2 as well as 
4 in the five-site case all possess this refiection symmetry 
(while the rest are anti-symmetric). This symmetry is 
not broken in the presence of nonlinear coupling between 
the sites. To achieve their dynamical generation we make 
use of a second property of the linear modes: the phase 
oscillation is half a cycle out of phase with the amplitude 
oscillation, i.e. the phase is at the extremum of an os- 
cillation, while the amplitude is at the mid-point. This 
means that with an appropriate choice of phase, the os- 
cillations may be initiated from the stationary state. To 
simulate the possible experimental conditions we apply a 
fixed phase shift between maxima of the periodic poten- 
tial, but extending to the boundaries of our simulation 
domain in the dimension perpendicular to the oscilla- 
tor chain. The phase perturbations for the three-site, 
four-site and five-site modes we excite are respectively: 
(0, e, 0), (0, e, e, 0) and (0, -e/2, e/2, -e/2, 0). The oscilla- 
tions of the particle numbers in each well, resulting from 
these phase perturbations, are shown in Fig. [3] for the 
three-site, four-site and five-site modes mentioned above. 
Note, Ni = Jy I'^pdr, where Vi signifies the spatial do- 
main of the i-th occupied well of the periodic potential 
V{x^ y). We see that in each case the oscillations are large 
in amplitude and regular, illustrating the persistence of 
the linear modes into the nonlinear regime. We show 
in Fig. [3fb,d,f) the dependence of the oscillation fre- 
quency on the phase perturbation e (and corresponding 
amplitude of oscillation), and see that as the amplitude 
increases the oscillation frequency begins to deviate from 
the linear oscillation frequency, as expected for a non- 
linear system. The spectrum also shows the appearance 
of a frequency subharmonic in the five site case, which 
increases with the perturbation, ultimately leading to a 
strong perturbation of the linear mode. 

As may perhaps be expected, not all of the linear 
modes are robust in the highly nonlinear regime. We 
focus here on some of the key instability/persistence fail- 
ure examples which may be observed in this coherent 
wave oscillator system. The first illustrated case exam- 
ple where such a mode does not persist has no paral- 
lel in analogue mechanical systems, and stems purely 
from the extended nature of the underlying system. This 
nonlinearity-induced 'cut-off' in possible mode excitation 




0.2 ^ 0.4 



FIG. 3: (color online). Mass/power oscillations per site for 
(a) three-site state, e = 0.27r (thick solid line: central site); 
(c) four-site state, e = O.Stt (solid lines: central sites); (e) 
five-site state, e = 0.27r (thick black line: central site; thick 
dashed lines: adjoining sites; thin dashed lines: outer sites). 
(b,d,f) Frequency dependence of resultant oscillations associ- 
ated with perturbation e for three-, four- and five-site cases, 
showing deviation from linear predictions (dashed lines) at 
large amplitude. Data presented is the Fourier transform of 
trajectories propagating to t = 200. 



occurs for all modes, however we consider the simplest 
two occupied site case. In a two-mode reduction, the 
Josephson-like oscillation between the two sites always 
exists for a small enough perturbation from the ground 
state configuration (TYq, A^o)- As can be seen in Fig. [2](a) 
a significant effect of the extended periodic lattice is that 
this oscillation ceases to exist, even in the infinitesimal 
limit, for stationary states with large enough chemical 
potential /i. Thus even though the stationary two-site 
state is stable, no localized internal mode of amplitude 
oscillations is possible in this regime of large ji. As seen in 
Fig. [3Ka) an appropriate phase perturbation (e, 0) does 
not lead to any coherent tunnelling dynamics between 
the sites. In this regime, it is evident that the simple 
two-site approximation breaks down, although the lin- 
earization captures this feature. 

The second considered pathway is commonly encoun- 
tered in mechanical systems, and consists of mode in- 
stability due to nonlinearity [2|, |3[. As an example we 
consider the first mode of the three-site case, which is 



analogous to the in-phase mode of two coupled pen- 
dula. We find this mode to be unstable for our coher- 
ent wave oscillator case, with the dynamics (see Fig. 
ID^b)) apparently pointing to an oscillatory instability. 
This can be attributed to a parametric instability due to 
the nonlinearity-induced coupling between the different 
modes in a way similar to what happens in the mechani- 
cal system analyzed in [3[. Here, linearization still offers 
the guiding principles of the relevant modes, but nonlin- 
earity destabilizes (some of) them, via mode coupling. 

The final pathway that we consider for the drastic 
modification of linear dynamics under strong nonlinear- 
ity is the breakdown of regular oscillations at large am- 
plitudes (i.e., the breakdown of linearized analysis alto- 
gether as a guiding principle for the dynamics). As we 
can see in Fig. lU^c), at large amplitude (corresponding 
to initial phase perturbation e = OAtt) the excitation 
of mode 4 in the five-site state ceases to be regular, as 
we also saw in the spectral profile of Fig. [3ff). Oscilla- 
tions persist, however the initial regular oscillation breaks 
up and some energy is lost through tunnelling to other 
wells. As an interesting aside, the reflection symmetry 
is still maintained throughout the irregular dynamics. If 
we look back to our discrete equations we see that for 
dynamics satisfying this reflection symmetry, chaos is in- 
deed only possible for five sites or more. The four site (re- 
flection symmetric) case amounts to effectively two sites 
but with two conservation laws, so no chaos is possible. 

Finally, we show that the oscillator framework we have 
introduced may be easily extended to genuinely higher 
dimensional lattice configurations. As an example we 
consider the simplest genuinely two-dimensional configu- 
ration, the four-site square shown in Fig. [5fa). Phase os- 
cillations of an initial vortex phase have been previously 
predicted in this configuration [23| (see also [2^ for the 
excitations of the lattice contour). Here, we illustrate 
the amplitude oscillations that may also be observed in 
this geometry. Following the symmetry arguments of the 
chain configuration we note that this system is symmet- 
ric when reflected through a line running through the 
centre and two sites, and also through a line running 
through the centre and bisecting the occupied sites. We 
find modes satisfying one of these symmetries are stable 
and we consider here the excitation of the diagonal mode. 
By imposing a phase shift on the two diagonally oppo- 
site sites (see Fig. [5fc)), we initiate a coherent tunnelling 
oscillation between the sites (see Fig. EJe)). These os- 
cillations correspond to a diagonal oscillation of phase 
and amplitude (Fig. [5fb,d)) and persist even in the large 
amplitude limit. We have observed such features (e.g., 
long-lived periodic oscillations) in more extended square- 
like states, however a detailed discussion of the complex 
oscillation patterns possible therein is beyond the scope 
of this work and will be analyzed elsewhere. 

In conclusion, we have presented some prototypical fea- 
tures within the rich potential of a framework for un- 
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FIG. 4: (color online). Examples of instabilities: (a) Absence 
of Josephson-lik e oscillations for perturbations of the station- 
ary state with /x = 7 (e = O.Itt); (b) Unstable first linear mode 
in three-site state (transfer of energy to the central site with 
e = iO.OSvr on outer sites); (c) Breakdown of mode 4 in the 
five-site case for e = OAtt. The lines follow the same labelling 
convention as in Fig. [H /i = 6 in (b,c). 



derstanding and predicting complex dynamics of wave 
systems in a periodic potential. We have shown that 
Josephson-like oscillations can still be observed, and 
thus do not require the spatial confinement imposed 
by a double-well potential. We have revealed that the 
Josephson-like dynamics may be understood as an oscil- 
lation mode of the stationary two-site state of the poten- 
tial. We have generalised this formalism to a chain of M 
coupled sites, revealing that such a chain has M — 1 linear 
oscillation modes. We have reduced the corresponding 
dynamics to a discrete dynamical system, showing that 
the observed oscillations may be understood as the nor- 
mal modes of this effective oscillator system. We have 
shown that some of the modes are robust even in the 
large-amplitude, nonlinear case, and we have provided a 
set of guidelines for their experimental excitation. On 
the other hand, we have also proposed a series of non- 
persistence/instability pathways that may arise, ranging 
from the disappearance of Josephson-like modes, to in- 
stabilities due to mode coupling and even to complex ir- 
regular dynamics beyond the mode description. We have 
shown that this formalism can be extended to higher di- 
mensions, and provides again a valuable theoretical tool 
for exploring oscillation modes therein. A systematic ex- 
ploration of such higher dimensional settings would be a 
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FIG. 5: (color online). Amplitude/phase oscillations for a 
two-dimensional four-site configuration, (a) Density and (c) 
associated phase with initial phase perturbation correspond- 
ing to e = 0.57r at t = (b,d) corresponding density and phase 
at t = 20 (values given by color-bars on right), (e) Associated 
power-per-site oscillation. 



natural step for future investigations. 
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